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ABSTRACT 

We have measured the Zeeman splitting of OH megamaser emission at 1667 MHz from five (ul- 
tra)luminous infrared galaxies ([U]LIRGs) using the 305 m Arecibo telescope and the 100 m Green 
Bank Telescope. Five of eight targeted galaxies show significant Zeeman-splitting detections, with 
14 individual masing components detected and line-of-sight magnetic field strengths ranging from 
~ 0.5 — 18 mG. The detected field strengths are similar to those measured in Galactic OH masers, 
suggesting that the local process of massive star formation occurs under similar conditions in (U)LIRGs 
and the Galaxy, in spite of the vastly different large-scale environments. Our measured field strengths 
are also similar to magnetic field strengths in (U)LIRGs inferred from synchrotron observations, imply- 
ing that milligauss magnetic fields likely pervade most phases of the interstellar medium in (U)LIRGs. 
These results provide a promising new tool for probing the astrophysics of distant galaxies. 
Subject headings: galaxies: magnetic fields — ISM: magnetic fields — magnetic fields — masers — 
polarization — radio lines: galaxies 



1. INTRODUCTION 

(Ultra)luminous infrared galaxies ([U]LIRGs) are a 
population of galaxies that emit far- infrared (FIR) radia- 
tion with energies comparable to those of the most lumi- 
nous quasars (logfXFiR/io) > H and 12 for LIRGs and 
ULIRGs, respectively). Nearly every ULIRG appears to 
have undergone a merger/interaction and contains mas- 
sive star formation and/or an active galact ic nucleus 
(AGN) induced by gravitational interactions. |Lc] (|2005f ) 
details very long baseline interferometry (VLBI) obser- 
vations of the 1667 MHz hydroxyl (OH) transition in the 
nuclear regions in (U) LIRGs that have revealed multiple 
masing regions with 1 < log(LoH/^o) < 4; these regions 
are known as OH megamasers (OHMs). Each OHM has 
a spectral linewidth of between 50 and 150 km s _1 ; when 
viewed by a single dish, these spectral components are 
superimposed. The 1667 MHz OHM flux density is al- 
ways a few to many times that of the 1665 MHz tran- 
sition, and in many cases t he 1665 MHz line is absent 
([Darling fc GiovanelU 12002): this is an interesting con- 
trast to the case of OH masers in the Gala xy in which the 
1665 MHz transition is usually dominant (jReid fc Moranl 
119881) . The star bursts and AGNs in ULIRGs create 
strong FIR dust emission as well as a strong radio con- 
tinuum; the OHMs are generally believed to be pum ped 
by the FIR radiation field (e.g., iRandell et aljfl995h al- 
thoug h collisional excitatio n may be important as wel l 
(e.g.. lLonsdale et all Il998l ). lLockett fc Elitzurl (|2008l ) 
have recently suggested that the 53 /^m OH pump lines 
in addition to line overlap of large (> 10 km s _1 ) turbu- 
lent linewidths can account for the observed dominance 
of the 1667 MHz transition in OHMs. They further argue 
that pumping due to FIR radiation can explain all ob- 
served main-line OH masers, both those in Galactic star- 
forming regions and those in OHM galaxies. Given the 
conditions that exist in ULIRGs and considering that so 
many OH masers in our Galaxy are associ ated with mas- 
sive star-forming regions ([Fish et al.ll2003l ). it is therefore 
not surprising that the entire OHM sample finds homes 



in LIRGs , strongly favoring the most FIR-luminous, the 
ULIRGs ((Darling fc Giovanellil[2002] ). 

The high gas and energy densities in ULIRGs make 
them natural locations to expect very strong magnetic 
fields. Much of the radio emission in ULIRGs is re- 
solved on scales of ~ 100 pc with VLA observations 
(|Condon et al.1 Il991h. High -resolution observations of 



Arp 220 (jRovilos et al.l f2003) show that the OHMs arise 
in this region as well. With this size scale and the ob- 
served radio flux densities, minimum energy arguments 
suggest volume averaged field strengths of ~ 1 mG (e.g., 
ICondon et all 119911 : [Thompson et all l2006| ). which are 
significantly larger than the ~ 10 fiG fields in normal 
spirals. The field strengths in ULIRGs cannot be much 
below a mG or else inverse Compton cooling would dom- 
inate over synchrotron cooling, making it energetically 
difficult to explain the radio flux densities from ULIRGs 
and the fact that ULIRGs lie on the FIR-radio correla- 
tion. The field strengths could, however, in principle be 
larger than the minimimum energy estimate if, as in our 
Galaxy, the magnetic energy density i s in approximate 
equip artition with the total pressure (jThompson et al.l 
2006J). The latter can be estimated from the observed 
surface density. CO observations of Arp 220 and several 
other systems reveal ~ 10 9 M P > of molecular ga s in the 
central ~ 100 pc (e.g.. iDownes fc Solomonlll998h imply- 
ing gas surface densities S ~ 1 — 10 g cm~ 2 , 10 3 — 10 4 
times larger than in the Milky Way (MW). The equipar- 
tition field scales as B cx S, implying that the mean field 
in ULIRGs could approach ~ 10 mG. 

Motivated by the above considerations, we carried out 
a survey of 8 (U) LIRGs searching for Zeeman splitting 
in OHMs. This paper presents our results, which repre- 
sent the first detections of extragalactic Zeeman splitting 
from an emission line and the first extragalactic detec- 
tions within an external galaxy proper. T he only pre- 
vious extragalactic dete ction was made by iKazes et al.l 
(|1991h and confirmed bv lSarma et al.l (|2005f) via absorp- 
tion of 21 cm emission in a high- velocity system towards 
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NGC 1275 (Per A). §[2]outlines what is known about each 
of our targets. In § [3] we describe the observations. In 
§U]we discuss the data reduction and calibration method. 
§§ [5] and [6] present a summary and discussion of the re- 
sults, respectively. 

2. SOURCE SELECTION 

In selecting our sam ple of targets from the comp i lation 
of all known OHMs bv lDarling k Giovanellil (|200Ct 120011 . 
2002), we chose the three simplest criteria possible. We 
selected 12 (U)LIRGs: (a) with the largest OHM peak 
flux densities; (6) whose discoverers did not regard the 
OHM detection validity as suspicious; (c) that are ob- 
servable from Arecibo, Puerto Rico, or Green Bank, West 
Virginia. Our sample incl udes 2 of only 3 known OH 
gigamasers {Loh > 10 4 L g : iDarling k GiovaneIlill2002f K 
Here we summarize what is known about each source and 
its OHM emission. 1 

IRAS F01417+1651: This LIRG is most commonly 
known 2 as III Zw 35 and has an optical heliocentric red- 
shift of z = 0.0274. It is a doub le galaxy system and 
is clas sified as a Seyfert 2 galaxy. IStavelev-Smith et al.1 
(|1987| ) present a single-dish spectrum from the Jodrell 
Bank MklA 76 m telescope showing emission from the 
1667 MHz transition at a velocity of 8262 km s -1 with a 
peak flux density of 240 mJy and a total velocity extent 
of 270 km s _1 at the 10% flux density level. The line 
profile can easily be seen to have at least three compo- 
nents. The 1665 MHz line is also weakly detected (~ 25 
mJy) and completely separated from the 1667 MHz emis- 
sion, with an estimated hyperfine line ratio (defined as 
Rh = f j 1667 dv/ J /i665 du, where the integrals repre- 
s ent the total flux den sity of each transition) of Rh ~ 9. 

iKilleen et all (|1996l ) observed III Zw 35 using the Aus- 
tralia Telescope Compact Array (ATCA). Their goal was 
an attempted detection of Zeeman splitting in the OHM 
emission. The peak flux density was 247.6 mJy and their 
sensitivity was 3.6 mJy. They observed no Zeeman split- 
ting and their model-dependent estimate for a 3a upper 
li mit on the line-o f-sight magnetic field was 4.0 mG. 

iDiamond et al.l (|1999f ) present global VLBI observa- 
tions of the OHM emission in III Zw 35. They label 
two regions of 1667 MHz OHM emission in the south (SI 
and S2) and three in the north (N1-N3), each region cov- 
ering about 20 mas and separated by 9 mas. They re- 
cover 60% of the single-dish flux density. Pihlstro m et al.l 
(|2001h performed simultaneous high-resolution observa- 
tions of the OHM emission in III Zw 35 using both the 
European VLBI Network (EVN; baselines between 198- 
2280 km) and the Multi-Element Radio-Linked Interfer- 
ometer Network (MERLIN; operated by Jodrell Bank 
Observatory with baselines between 6.2-217 km). The 
map of the 1667 MHz emission shows two compact re- 
gi ons coincident w i th the northern and southern sources 
of IDiamond et al.1 (|1999h connected by two bridges of 

1 Unfortunately, we only obtained usable data for 2 of the 6 
sources observed at Green Bank, therefore we only provide source 
descriptions for the 8 (U)LIRGs for which we have presentable 
results. 

2 As a shorthand, when referring to source names in the text by 
their IRAS designation, we shall henceforth use the right ascension 
designator only; we shall refer to 01417, 10038, and 15327 by their 
more common designators, III Zw 35, IC 2545, and Arp 220, re- 
spectively. We shall retain the full IRAS designation in figure and 
table captions as well as section headings. 



weaker, more diffuse emission. In total, 80% of the single- 
dish flux density was recovered. A velocity gradient of 
~ 1.5 km s _1 pc _1 is observed from the southern to the 
northern regions and is evident in the diffuse component. 
The emission is modelled as a torus of multiple mascr 
clouds inclined at 60° ; the compact OHM emission would 
be seen at the tangent points where a few clouds could 
be superimposed in such a fashion that strong OHM 
emission would be produced from the foreground clouds 
amplifying those in the background. At the front and 
back of the torus, the emission would be weak because 
the path lengths through the to rus are small and the 
clouds are less li kely to overlap ([Pihlstrom et al.1 [20011 ; 
iParra et al.ll2005h . 

IRAS F10038-3338: Also known as IC 2545, this 
LIRG is a set of interacting galaxies at z = 0.0341. A 
single-dish spectrum made using the Park e s 64 m tele- 
scope is presented bv lStavelev-Smith et al.l (|1992T ) show- 
ing 1667 MHz OHM emission centered at 10093 km s" 1 
with a full width at half maximum (FWHM) velocity 
range of 63 km s _1 and a peak flux density of 315 mJy. 
Likely due to its low declination, there have been no 
V LBI observations of this source despite its brightness. 

Kille en et alj (|1996f ) present an ATCA spectrum with 
much better sensitivity (5. 4 mJy in Stokes I) an d veloc - 
ity resolution than that of IStavelev-Smith et al.l (l992); 
the 1667 MHz emission contains five narrow peaks super- 
imposed on a broad emission component. The brightest 
component has a peak fl ux density of 260 m Jy and a 
velocity of 10097 km s _1 . IKilleen et all (|1996h failed to 
detect Zeeman splitting and estimated that the line-of- 
sight field should be less than 4.3 mG. 
IRAS F10173+0829: The only single-dish observa- 
tions of the OHM line emission in this LIRG at z = 
0.0480 were m ade with the 305 m Arecib o telescope and 
are detailed in lMirabel k Sanders! (|1987h . There are two 
distinct peaks in the profile of the 1667 MHz emission 
with a separation of about 100 km s _1 , the dominant 
peak having a velocity of 14720 km s _1 with a FWHM 
of 39 km s _1 and a peak flux density of about 105 mJy. 
The 1665 and 1667 MHz lines are well separated with a 
hyperfine line ratio Rh = 14.6. 

MERLIN observations made bvlYul (|2004l I2005D show 
roughly 50 maser spots distributed into three clumps, la- 
belled east, central, and west, over an area of 1'A x Of 6. 
The spots within each clump are distributed along a line, 
with each of the three lines having a different direction; 
E3 (2005) proposes (without much justification) that the 
spots may be distributed along a warped circumnuclear 
torus seen edge-on. The OHM emission is seen only at 
the 1667 MHz transition and is coincident with the in- 
frared central position. 

IRAS F11506-3851: Also known as ESO 320-30, 
this LIRG is classified as an H II galaxy at z = 
0.0108. A Parkes sin gle-dis h spectrum is presented by 
IStavelev-Smith et all (fl992T ) showing 1667 MHz OHM 
emission centered at 3103 km s^ 1 with a FWHM ve- 
locity extent of 87 km s _1 and a peak flux density of 
105 mJy. There is neither enough sensitivity nor band- 
width to clearly discern any 1665 MHz emission. There 
are no interferometric observations of this source: like 
IC 2545, the low declination of 11506 would hinder any 
attempted VLBI observations. 

IRAS F12032+1707: A gigamaser discovered at 
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Arecibo by iDarling fe GiovaneliU (|2001h . The host ob- 
ject, a ULIRG at z = 0.2170, has been classified as 
a LINER-type AGN (|Veilleux et al.|[l999h . The OHM 
emission spans almost 2000 km s _1 with a redshifted 
high-velocity tail and a mean flux density of roughly 9 
mJy. The 1665 and 1667 MHz lines are impossible to dis- 
tinguish and clearly blended. A very narrow and bright 
component is seen at 64500 km s _1 with a peak flux 

d ensity of 16.3 mJ y. 

Pihl strom et al.l (|2005l ) used the Very Long Baseline 
Array (VLBA) to show that the OHM emission is con- 
fined to an area of 25 x 25 mas. All the single-dish 
flux density was recovered. They were able to clearly 
identify five peaks in their Stokes / spectrum that cor- 
responded with Darling's single-dish spectrum. By aver- 
aging channels around each peak, they found the mascr 
components were spatially separated and aligned roughly 
north-south, implying an ordered velocity gradient. No 
continuum emission was detected, implying that the con- 
tinuum emission is resolved out on scales less than 75 
mas. 

IRAS F12112+0305: This ULIRG is classified as a 
LINER-type AGN and is an interacting pair of galax- 
ies at z — 0.0730. The only information concerning the 
O HM emission in t his ULIRG is listed in tabular form 
in IStavelev-Smith et al.l (|1992l ): no spectrum has been 
published. The 1667 MHz line was measured at a veloc- 
ity of 5540 km s _1 and no information about the 1665 
MHz tran sition is published. It s 1667 MHz flux density 
is listed in IDarling fc GiovaneliU (|2002l ) as 45 mJy. There 
are no VLBI observations of the 1667 MHz OHM emis- 
sion for this source. 

IRAS F14070+0525: Discovered by iBaan et all 

(1992), this gigamaser is the most distant OH M at a 
redshift of z = 0.2644. IDarling fc Giovanellil ((20021 ) 
redetected this source in their survey. The OH lines 
are so wide and blended (1580 km s _1 at 10% peak 
flux density) that it is impossible to identify any 1665 
MHz emission. The spect ral line profile measured by 
IDarling fc GiovaneT li (2002) has a peak flux density of 8.4 
mJy and sho ws no signifi c ant ch an ges since the original 



detect ion by IBaan et all (|1992t ). IDarling fc Giovanellil 



(2002) suggest that many peaks in the profile are likely 
the result of many masing nuclei within the host ULIRG, 

which is classified as a Seyfert 2. 

VLBA observations made by IPihlstro m et al.l ([2005) 
recovered less than 10% of the single-dish flux density 
and only two of the many single-dish spectral peaks were 
detected. Most of the single-dish emission is therefore 
diffuse. The spatial extent of the VLBA emission is con- 
fined to 10 x 10 mas. 

IRAS F15327+2340: This is perhaps the most well- 
know n ULIRG and is b etter known as Arp 220 or IC 
4553. IBaan et al.1 (| 1982ft discovered the OHM emission 
using Arecibo and list the single-dish properties of the 
OHM emission as having a velocity of 5375 km s _1 and 
FWHM velocity extent of 108 km s _1 . The spectrum 
clearly shows that the 1665 and 1667 MHz transitions 
are distinct with a hypcrfinc line ratio Rh — 4.2 and a 
p eak flux dens i ty of 320 mJy. 

iSmith et al.l (|1998ft used global VLBI continuum imag- 
ing at 18 cm to show that the high brightness tem- 
perature core of Arp 220 is composed of multiple com- 
pact sources, which they interpret as luminous radio 



supernovae (RSNe). These RSNe are not coincident 
with the comp act 1667 MHz OHM spot s discovered by 
lLonsdale et al.l (|1998l ). More recently, lLonsdale et al.1 
(2006) have used high-sensitivity 18 cm VLBI observa- 
tions of the nuclei to detect 4 previously unseen sources 
in a one year perio d, supporting the RSNe interpretation. 
iParra et al.l (|2007| ) have made the first multiwavelength 
observations of these compact sources; they identify a 
fr action of these sourc es to be supernova remnants. 

iRovilos et all (|2003D present MERLIN maps of the 
18 cm continuum and OHM line emission in Arp 220; two 
components are seen roughly 1" apart, e ach coinciding 
with a nucleus imaged in the infrared by iGraham et al.l 
(1990). The OHM emission is resolved into one com- 
ponent aligned with the eastern continuum feature and 
two components that are ali gned north to south s trad- 
dling the wes t ern n ucleus. lLonsdale et al.1 (|1998l ) and 
IRovilos et al.l (|2003l ) present global VLBI spectral line 
maps that show that the OHM emission is resolved 
into multiple compact spots. The northernmost features 
in both the eastern and western nuclei form elongated 
ridges. 

3. OBSERVATIONS 

In 2006 February we used the L-band wide receiver of 
the 305 m Arecibo 3 telescope in full-Stokes mode in an 
attempt to detect Zeeman splitting of the 1667 MHz OH 
transition in the 6 positive-declination sources listed in 

§12 

Since the spatial extent of each source is much smaller 
than Arecibo's 3.'3 beam, our observing method was to 
simply spend equal time at on-source and off-source po- 
sitions. In this position-switching scheme, we alternated 
between 4 minutes on source and 4 minutes at a reference 
position having the same declination as the source and 
a right ascension 4 minutes east of the source. In this 
way, the hour angle ranges of the source and reference 
observations were nearly identical. Our integration time 
was one second, allowing us to remove short-term radio- 
frequency interference (RFI). The total integration time 
for each source is as follows: 5.4 hours for III Zw 35, 2.7 
hours each for 10173 and 12032, 3.1 hours for 12112, 4.6 
hours for 14070, and 5.9 hours for Arp 220. We config- 
ured the correlator to produce 4 spectra per integration: 
one 6.25 MHz bandpass centered on the mean of the 1665 
and 1667 MHz transitions; 2 narrow bandwidths (either 
3.125 or 12.5 MHz, depending on the velocity extent of 
the source) centered on the 1665 and 1667 MHz transi- 
tions, respectively; and one wide bandwidth (either 12.5 
or 25 MHz) centered on the mean of the 1665 and 1667 
MHz transitions. We calibrated the Mueller matrix for 
Stokes parameters using the standard Arec ibo technique 
(jHeiles et al.l 120011 : iHeiles fc Trorandll2004f ) of observing 
spider scans on the linearly-polarized continuum sources 
3C 138 and 3C 286. 

In 2005 December we used the L-band receiver of the 
100 m Robert C. Byrd Green Bank Telescope 4 (GBT) to 

3 The Arecibo Observatory is part of the National Astronomy 
and Ionosphere Center, which is operated by Cornell University 
under a cooperative agreement with the National Science Founda- 
tion. 

4 The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 
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observe an additional six OHM galaxies. All but the two 
sources in § [2] at negative declinations were affected by 
insidious RFI that left our data corrupted beyond salva- 
tion (For completeness, the sources that were obliterated 
by RFI were IRAS F12540+5708, IRAS F13428+5608, 
IRAS F17207-0014, and IRAS F20100-4156.) We used 
two observing methods for each source: position switch- 
ing (as described above) and least-squares frequency 
switching (LSFS; see lHeilesll2007l for the details of this ob- 
serving method and its corresponding reduction scheme) . 
The LSFS method was used to accurately derive the gain 
for each integration; the data were then combined in 
the standard way using the off-source, position-switched 
spectra. Our off-source positions were 23 minutes east 
of each on-source position in order to cause the GBT to 
track as closely as possible the path of our on-source ob- 
servations. We used a 12.5 MHz bandwidth and 9-level 
sampling for all observations. 

We accumulated 4.0 hours of RFI-free, on-source inte- 
gration time for IC 2545 and 5.8 hours for 11506. As we 
did at Arecibo, we observed 3C 286 using spider scans 
in order to calibrate the L-band Mueller matrix at the 
GBT. 

4. DATA REDUCTION 

The complex Stokes / (see footnote 5 ) line shape in each 
of the maser sources is a composite of many narrow maser 
lines at various velocities spread about the systemic ve- 
locity of the system. Therefore, we chose to least-squares 
fit each line profile with multiple Gaussian components. 
Without VLBI observations, it is impossible to attribute 
any particular velocity or width to a Gaussian compo- 
nent within the profile. The only method available to us 
for assessing a possible field strength from each Stokes V 
spectrum was to decompose each I profile into the fewest 
number of Gaussian components that would yield rea- 
sonable residuals while also allowing enough components 
to reproduce the multiple splittings in the V spectrum. 
None of the parameters in our multiple-component Gaus- 
sian fits were held fixed. We discuss our method in more 
detail in § |4~! 

4.1. Calibration 

The derived Mueller matrix was applied to all OHM 
observations to correct the polarization products and ob- 
tain the pure Stokes spectra for each observed source. We 
converted from antenna temperature to flux density by 
assuming the antenna gain to be 10.0 K Jy _1 at Arecibo 
and 2.0 K Jy" 1 at the GBT; these gains were estimated 
from observations of standard flux density calibrators. 

We follow the IAU definition for Stokes V, namely 
V = RHCP - LHCP, where RHCP is the IEEE def- 
inition of right-hand circular polarization. 6 We deter- 
mined the sense of Stokes V at Arecibo by observing 
the highly circularly polarized 1665 MHz Galactic maser 
W49(OH); the result i s con sistent with the measure ments 
of iRoeers et "all (|1967ft and lColes & Rumsevl (|1970t ). The 

5 We use the classical definition of Stokes I, which is the sum 
(not the average) of two orthogonal polarizations. Thus, stated 
Stokes / flux densities are twice those listed in § [2] and in other 
catalogs. 

6 Defined as a clockwise rotation of the electric vector along the 
direction of propagation. 



sense of Stokes V has not yet been determined for the 
GBT Autocorrelation Spectrometer, which only began 
functioning with full-Stokes capability months prior to 
these observations. 

4.2. RFI Removal 

We examined each set of spectra, both off source and 
on, for RFI and rejected suspicious-looking data, which 
constituted only a few percent for only two sources, III 
Zw 35 and Arp 220. The other sources were completely 
free of RFI except for occasional monochromatic signals 
whose topocentric frequencies are constant; fortunately, 
most of these fall off the OHM lines. 

For one source, 10173, the monochromatic RFI fell on 
the OHM line. We observed this source over several days, 
during which the changing Doppler shift moved the RFI 
across part of the OH spectrum. For each day the RFI 
was a sharp spike with the usual ringing sidelobes. We 
Hanning smoothed each spectrum, which eliminated the 
ringing, and interpolated across each day's spike, which 
effectively removed the RFI. 

The source III Zw 35, whose OH lines are centered 
near 1622.5 MHz, was highly contaminated by RFI that 
probably arises from the Iridium communications satel- 
lites. The RFI consists of a spikey pattern that repeats 
periodically across the spectrum at about a 0.33 MHz 
interval. It was impossible to obtain reasonable results 
by averaging data. However, by taking medians instead 
of averaging, the spectra look quite good and the RFI 
is reduced to levels of about 20 mJy in Stokes I and 3 
mJy in Stokes V — levels that are considerably smaller 
(< 1%) than the OHM spectral features. 

4.3. Bandpass and Gain Correction 

We correct our spectra for the intermediate-frequency 
bandpass. Both the Arecibo L-band wide receiver and 
the GBT L-band receiver are dual-polarized feeds with 
native linear polarization. For Stokes I and Q, we divide 
each of the two linear polarization spectra (XX and YY) 
by its associated off-source spectrum; then we add the 
results to obtain Stokes I and subtract them to obtain 
Q. To generate Stokes U and V, we combine the cross- 
correlation spectra (XY and YX) having divided by the 
square root of the product of the off-source XX and YY 
spectra. 

We always show difference spectra: on-source minus 
off-source. Normally, on- and off-source spectra are com- 
bined by subtracting the latter from the former. If the 
two spectra have equal noise a, then the noise in the 
difference is y/2a. Our off-source spectra have no fine- 
scale frequency structure, so we can reduce the noise by 
smoothing. We use a Fourier technique to smooth the off- 
source average spectrum. By zeroing lags at high delays 
in the autocorrelation function of the average off-source 
spectrum and then Fourier transforming, we nearly elim- 
inate the noise contribution from the off-source spectrum 
while retaining the shape of the bandpass. This reduc- 
tion in noise is particularly important for the polarized 
Stokes parameters, which are weak. 

As mentioned above, Stokes U and V are obtained via 
the cross-correlation products XY and YX: this insu- 
lates them from system gain fluctuations. However, since 
Stokes Q is the difference between the two native linear 
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polarizations, it is susceptible to gain fluctuations. We 
defer the discussion of linear polarization to § 15.21 

After gain and bandpass correction, the on- and off- 
source spectra were averaged separately and combined 
to yield the final average Stokes spectra. 

4.4. Fitting Gaussian Components to Stokes I Profiles 

Fitting multiple Gaussians to complicated spectral 
profiles carries a significant degree of subjectivity because 
the fits are nonlinear. Generally, nonlinear fits require 
beginning from initial "guessed" parameters and letting 
the f it converge with successive iterations ([Press et al.l 
1992). Nonlinear fits usually have multiple x 2 minima 
and the particular minimum selected depends on the 
initial guesses, which in turn depend on the subjective 
judgement of the person doing the fitting. Therefore, we 
outline the following guidelines that were used for select- 
ing the initial guesses for Gaussian components in the 
fits to the Stokes / spectra: 

1. For each peak (i.e., local maximum) in /, we in- 
cluded a single Gaussian component whose three 
parameters (flux density, central frequency, and fre- 
quency width) were visually estimated. 

2. Many / peaks are distinctly asymmetric. We fitted 
these asymmetries by including one or two Gaus- 
sian components with visually-estimated parame- 
ters in addition to the central component of guide- 
line (1). 

3. The Gausssian components estimated in guidelines 
(1) and (2) are usually fairly narrow and lie on top 
of one or two underlying broader lines — core-halo 
structure. We included one or two broad Gaussian 
components to represent these broader lines. 

4. With all of the above, our goal was to use the fewest 
number of Gaussian components that would yield 
reasonable residuals. 

It was straightfoward to apply the above guidelines to 
the sources III Zw 35, IC 2545, 11506, 12112, 14070, and 
even Arp 220, for which we fitted 18 components. The 
sources 10173 and 12032 are somewhat more complex. 
In §§ 15.1.1^5X51 we describe how we applied the above 
selection guidelines when appropriate. 

We stress that our Gaussian-component representa- 
tions are not unique. Particular problems include: 

1. In guideline (1) above, noise prevents us from iden- 
tifying weak components. This introduces a sensi- 
tivity cutoff. Noise also prevents us from distin- 
guishing two or more closely-spaced blended real 
components from a single broader component. Be- 
cause we favor choices with the fewest number of 
components, this introduces a bias towards wider 
components. 

2. In guidelines (2) and (3) above, whether to rep- 
resent a peak needing multiple Gaussians by an 
asymmetry or core-halo structure can be extremely 
subjective. For example, the combination of two 
narrow Gaussians separated by a fraction of their 
FWHMs can closely mimic the combination of a 
broad and a narrow Gaussian with roughly the 
same centers. 



In summary, for most sources our Gaussian fits fol- 
low our fitting guidelines in a reasonably straightforward 
fashion; if the fits were done by other people who followed 
these guidelines, the components would be mostly repro- 
duced. However, Gaussian component fitting has uncer- 
tainties as mentioned above, particularly when compo- 
nents are blended and signal-to-noise (S/N) is low. 

5. RESULTS 

For each source, the Stokes / spectrum exhibits a fairly 
broad, relatively smooth underlying profile for the OHM 
emission on top of which small bumps from individual 
masers can be seen. VLBI studies show that the un- 
derlying profile arises from spatially extended OH emis- 
sion and, sometimes, an assembl y of many masers that 
are not individually recognizable (Pihls trom et al.l 120051 : 
iDiamond et al.llT999l ). We fit the Stokes / spectrum for 
each source with a series of Gaussian profiles. For Arp 
220, 12112, and 12032 we also fit a first-degree polyno- 
mial (12032 required a second-degree polynomial in ad- 
dition), since these profiles exhibit broad wings. A few 
sources have a large number of discernable individual 
masers: for example, we used 18 Gaussians components 
for Arp 220. 

We examined circular polarization for each source, and 
linear polarization for the six Arecibo sources only. Five 
of the eight sources exhibit significant circular polariza- 
tion that is interpretable as Zeeman splitting, particu- 
larly for the recognizable individual maser components. 
For four of the sources, there is evidence that the mag- 
netic field reverses direction between OHM spots within 
the source. 

We see detectable linear maser polarization in two 
sources (possibly four) and are able to estimate Faraday 
rotation in both. We present all spectra as a function 
of frequency as viewed in the heliocentric frame. Since 
all OHMs are extragalactic sources, OHM spectra are al- 
most always presented versus optical heliocentric velocity 
t>0, which is conventionally defined as 



c 



Hi 
v 



1 



Zq 



(1) 



where c is the speed of light, vq is the rest frequency 
(which is taken to be 1667.359 MHz for OHMs since this 
transition always dominates the 1665.4018 MHz transi- 
tion), v is the observed frequency, and Zq is the redshift 
of the maser. 

First, in § 15.11 we present the circular polarization re- 
sults for each source in addition to describing the total 
intensity properties. We present the linear polarization 
results for each source observed with Arecibo in § 15.21 
For sources in which the 1665 and 1667 MHz emission 
are separable, we calculate the hyperfine ratio Rr- 

5.1. Circular Polarization and Line-of-Sight Magnetic 

Fields 

For the usual case in which the Zeeman splitting is 
small compared to the linewidth, the Stokes V spectrum 
is given by 

(2) 



V=[^\( d S\hB 



dv J 



where B\\ is the line-of-sight component of the magnetic 
field at the OHM and b is known as the splitting coeffi- 
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cient 7 (|Heiles et alj|1993l ). equal to 1.96 Hz fiG^ 1 for the 
OH 1667.359 MHz transition; 8 the factor v/vq, equiva- 
lent to (1 + z Q ) _1 , accounts for the frequency compres- 
sion of redshifted lines. In order to derive a magnetic 
field strength, we need to least-squares fit the Stokes V 
spectrum with the functional form of equation (|2|). As 
is the custom in radio Zeeman work, we add a term on 
the right that is linear in Stokes I to account for leakage 
of / into the measured V. For Bu > (by convention a 
positive magnetic field points away from the observer), if 
Stokes V is plotted as a function of frequency, V will be 
positive on the low-frequency side of a Stokes I emission 
line. 

We solved equation ([2]) in two ways. In one, we simul- 
taneously fitted Stokes V for multiple Gaussian compo- 
nents (selected as outlined in § I4.4|) to derive separate, 
independent magnetic fields for each Gaussian. In the 
other, we chose a limited range in frequency is, either 0.1 
or 0.25 MHz, and fitted for Bu for the center of this range, 
positioning the center sequentially at each spectral chan- 
nel to obtain Bu as a function of frequency; we refer to 
this as the B(v) fit. The former method is appropriate 
for individual masers, while the latter is more suitable for 
the broad component. We plot the B(v) results only for 
sources for which the results provide additional insight. 

Performing these fits requires the calculation of dl/dv. 
The I profiles are somewhat noisy and the frequency 
derivative is often very noisy. This means that traditional 
least-squares fitting cannot be used, because it a s sumes 
no error in the independent variables. ISault et alj (|1990h 
discuss this and suggest using a generalized maximum 
likelihood technique. We choose the much simpler ap- 
proach of using our multiple-Gaussian fit to the Stokes I 
spectrum as the independent variable: it has no noise, so 
it satisfies the requirements of the conventional method 
of least squares. 

We present two vertically stacked plots for each source 
(Figures [THE]) . In the top panel, the Stokes / spectrum 
is plotted as a solid line and the profiles of the Gaus- 
sian components are plotted as dashed lines. The resid- 
uals (the difference between the data and the composite 
Gaussian fit) are plotted with enhanced vertical scale as 
a solid line near the middle height of the panel. Scale 
bars are plotted on both the residuals and the baseline 
of the spectrum to the right of the OHM emission; the 
height of each scale bar corresponds to the labelled flux 
density. The bottom axes of both plots show heliocen- 
tric frequency and the top axis of the top plot displays 
the optical heliocentric velocity. In the bottom panel, 
the Stokes V spectrum is plotted as a solid line and the 
dashed line represents the best fit to equation @. The 
integers located between the top and bottom panels label 
the number of each Gaussian component as assigned in 
the corresponding tabular summary and are positioned 
at the central frequencies of each component. 9 The dis- 

7 The splitting coefficient is directly proportional to the Landc 
g-factor for the transition: b = 2gfio/h, where fio is the Bohr 
magneton and h is Planck's constant. 

8 IModiaz et al.l 1120050 made a valiant effort to detect Zeeman 
splitting of 22.2 GHz H2O megamasers in NGC 4258 using the 
VLA and the GBT, but the splitting coefficient for this hypcrfinc 
transition is nearly 1000 times weaker than that of the 1667 MHz 
OH transition. 

9 Where multiple labels overlap, the font size has been reduced 



played spectra and residuals are smoothed over 7 chan- 
nels for every source. 

We assume that each Gaussian represents an emission 
component for the 1667 MHz transition. Arp 220 and 
12032 have line profiles that are too complex for the 1665 
and 1667 MHz lines to be distinguished. This introduces 
some uncertainty in our Zeeman splitting interpretations 
in §§ 15.1.51 & 15.1.81 For sources where the 1665 MHz 
transition is not blended with the 1667 MHz emission, 
we present the spectra showing both transitions in § 15.21 
and calculate the hyperfine line ratios. In all cases the 
1665 MHz transition was too weak for Zeeman splitting 
to be detected even if observed in the 1667 MHz line. 

5.1.1. IRAS F01417+ 1651 (III Zw 35) 

As we mentioned in § 14.21 our observations of III Zw 
35 suffered severe RFI that we were able to greatly re- 
duce by combining the data using medians instead of 
averaging. There remains a spikey pattern that repeats 
periodically across the spectrum at an interval of ~ 0.33 
MHz. Remarkably, this spikey pattern is restricted to an 
8 MHz wide interval centered almost exactly on the OHM 
lines. The spikey pattern appears in both the on-source 
and off-source spectra, so we regard this as terrestrial 
interference. Despite the RFI, Figure [T] shows that both 
/ and V are well-detected. In fitting Gaussians we are 
conservative because we realize that the RFI may have 
contaminated the line shape. In particular, the 0.33 MHz 
intervals happen to fall close to the two peaks in Stokes 
I. 

Table [T] lists the parameters of the Gaussian compo- 
nents that best fit the Stokes I spectrum. Column (1) 
lists the zero-based component number. Column (2) lists 
the peak flux density of each component in mJy and the 
corresponding uncertainty. Column (3) lists the central 
heliocentric frequency of each component in MHz and the 
corresponding uncertainty. Column (4) lists the FWHM 
of each component in MHz and the corresponding uncer- 
tainty. Column (5) lists the optical heliocentric velocity 
corresponding to the central frequency of each compo- 
nent. Column (6) lists the derived line-of-sight magnetic 
field in mG for each component and the corresponding 
uncertainty. 

Appying our guidelines from § 14.41 yielded the 5 Gaus- 
sian components shown in Figure [T] and listed in Table 
[TJ There are 2 distinct peaks in the OHM emission. The 
peak nearest 1622.8 MHz is asymmetric in such a way 
that 2 narrow components needed to be added near this 
peak in order to minimize the residuals. The overall pro- 
file has a core-halo structure, with the 4 narrow compo- 
nents lying on top of a broader component. 

The Stokes V spectrum (5 rms = 0.97 mJy) in Figure [1] 
shows prominent features that are fitted reasonably well 
by the 5 Gaussian components, with Zeeman splitting 
yielding significant fields in three Gaussians: Gaussian 
has Bu — 2.9±0.2 mG, and Gaussians 3 and 4 have fields 
of — 2.7 ± 0.1 and —3.6 ± 0.3 mG, respec tively. Thus the 
field r everses from one peak to the other. iPihlstrom et al.l 
(|2001h present 13 spectra from the 1667 MHz OHM emis- 
sion of III Zw 35 that they mapped using the EVN. These 
data show clearly that the 8215 and 8240 km s _1 compo- 

and the labels stacked corresponding to their associated flux den- 
sities. 
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nents (Gaussians 3 and 4) arise from the southern peak, 
while the 8312 km s _1 component (Gaussian 0) is as- 
sociated with the northern peak. This provides clear 
evidence that the reversal is arranged with the magnetic 
field pointing away from us in the north and towards us 
in the south. 

The panels on the right side of Figure [1] show results 
relevant to the B(v) fit. The top panel shows the compos- 
ite Gaussian-fitted (noise-free) Stokes I spectrum. The 
center panel shows the measured Stokes V spectrum as a 
solid line; the dashed line represents the Stokes V spec- 
trum that would be produced by a uniform line-of-sight 
magnetic field of 1 mG — this is obtained from equa- 
tion by setting Bu = 1 mG and using the derivative 
of the composite Gaussian shown in the top panel. The 
bottom panel shows the B{v) fit — the derived B\\ as a 
function of frequency — as described in § 15.11 There is 
a clear systematic pattern, with the field reversing sign 
from one peak to the other. The estimated field strengths 
a re also consistent w ith t he Gaussian fits. 

iParra et all (|2005h and lPihlstrom et all (|2001h present 
models of the OHM emission in III Zw 35 as a clumpy, 
rotating starburst ring at an inclination of 60°, with an 
inner radius of 22 pc and a radial thickness of 3 pc. 
Both the Gaussian and B(v) analyses suggest that an az- 
imuthal magnetic field is embedded within this starburst 
ring such that the field points towards us at the south- 
ernmost tangent point and away from us at the north- 
ernmost tangent point. IParra et all (|2005h estimate that 
the OHM clouds would be magnetically confined by a 
magnetic field of orde r ~ 10 mG. 

IKilleen et aT] (|1996l ) used the ATCA to observe the 
OHM emission in III Zw 35. Their sensivity of S Tma = 3.6 
mJy in Stokes V was not sufficient to detect the Zeeman 
splitting of the 1667 MHz line. 

5.1.2. IRAS F10038-3338 (IC 2545) 

As the residuals in Figure [2] show, the 1667 MHz OHM 
emission from IC 2545 is fitted extremely well by 5 nar- 
row Gaussian components and one broad one. Using our 
prescription from § 14.41 we see that there are 4 distinct 
peaks. The peak nearest 1613.1 MHz has an asymmetry 
that can be represented with a single extra narrow com- 
ponent. A broad component represents the evident core- 
halo structure of the OHM emission profile. It is unclear 
if the emission feature at 1611.7 MHz corresponds to 1665 
MHz emission or redshifted 1667 MHz emission. The 
Stokes I flux density and line profile have no t changed 
since the observations of IKilleen et all <fl996). Tabled 
lists the Gaussian fit parameters shown in Figure |U The 
Stokes V spectrum has an rms noise of S TmB = 0.7 mJy. 
There are three clear detections: Gaussian 1 probes a 
field of —1.8 ± 0.3 mG; Gaussian 2 is fitted by a field 
of — 11.3 ± 1.2 mG; and Gaussian 5 shows a reversal in 
sign with a field of 1.7± 0.3 mG. Since no VLBI observa- 
tions exist for this LIRG, nothing can be said about the 
structure of the field reversal. 

5.1.3. IRAS F10173+ 0829 

As mentioned in § 14.21 we used Hanning-smoothcd 
spectra when least-squares fitting this source because of 
RFI. We increased the derived unc e rtainties in Table CD 
by the appropriate factor of y/8/3 (|Killeen et alj 119961 
Table Al). 



We fit Stokes I with 7 Gaussians as shown in Figure 
[3] and listed in Table ED Using the selection guidelines of 
§ 14.41 we required 4 narrow components to sufficiently fit 
the extremely asymmetric peak near 1589.3 MHz. The 
fit to Stokes / yielded reasonable residuals by including 
a single extremely broad component. This source repre- 
sents a case where the profile structure is too complex 
to be modelled by our straightforward fitting guidelines: 
we regard the derived components as highly suspect re- 
gardless of the quality of the fit. 

Detectable portions of the Stokes V spectrum are re- 
stricted to the strongly peaked / line, where V oscillates 
rapidly. Due to the asymmetry in this profile, we found 
it impossible to obtain a set of Gaussians that gives good 
fits to both / and V. Our final fit reflects a compromise 
between large residuals and a larger number of Gaus- 
sians. While it is clear that Zeeman splitting has been 
detected in this source, we would need to add many more 
narrow, weak Gaussian components to our Stokes / fit 
in order to obtain statistically significant magnetic field 
derivations from the fit to Stokes V (which has an rms 
noise of S lms = 0.8 mJy). Without any VLBI observa- 
tions or any other physical motivation for adding such 
components, the best we can do is present the evidence 
for Zeeman splitting without estimating field strengths 
for individual maser components. 

Figure CO right panel, shows plots relevant to the B{v) 
fit. The dashed line in the middle panel shows the Stokes 
V spectrum that would be expected for a uniform B\\ = 2 
mG. The bottom panel shows the B(v) fit. There is 
a clear systematic pattern, with the field reversing sign 
from one side of the peak to the other. This reversal is 
not revealed by the Gaussian fits because there are no 
narrow Gaussians on either side of the peak and because 
the peak itself is represented by a single Gaussian (com- 
ponent 2). 

Since no sp ectr a l informatio n is presented in the MER- 
LIN maps of [Yu] (|2004l . 12005( 1 . it is impossible to asso- 
ciate any of our Gaussian components with OHM spots 
in 10173. 

5.1.4. IRAS F11506-3851 

The top panel of Figure [4] shows that the 1667 MHz 
Stokes I emission is fitted quite well by 6 narrow Gaus- 
sian components. The parameters for each component 
are listed in Table |H Gaussian components and 3 have 
derived magnetic fields that look significant when judged 
by their formal errors. However, given the quality of the 
Stokes V spectrum (S lms = 0.4 mJy), we have no confi- 
dence in either result. The 1665 MHz emission is clearly 
completely separated from the 1667 MHz emission, yield- 
ing a hypcrfinc ratio Rh = 4.9. 

5.1.5. IRAS F12032+ 1707 

The OHM emission from this gigamaser has an ex- 
tremely wide extent; it is impossible to distinguish the 
1665 and 1667 MHz emission. We fit its Stokes / spec- 
trum with 13 Gaussians as shown in Figure [5] with pa- 
rameters as listed in Tabled Due to the complexity of 
this profile, this is an extremely difficult case to apply our 
component selection guidelines to. We first identify 8 lo- 
cal maxima marked in Figure [5] by components 0, 2, 3, 5, 
6, 8, 10, and 12. The peaks near 8 and 10 display shoul- 
ders and are clearly blended with narrower components: 
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we added one additional component to each (components 
7 and 9, respectively). We represent a broad shoulder, 
not quite intense enough to produce a local maximum, 
near 1373 MHz by component 11. Finally, we add two 
broad halo components, numbers 1 and 4, to minimize 
the residuals. While the prescription of § 14. 41 allows us to 
select these 13 components somewhat straightforwardly 
— even producing respectable residuals — visual inspec- 
tion of Figure [5] suggests that our model simplifies and 
glosses over the innate complexity of this OHM profile. 
Without VLBI observations, there is little that can be 
done to improve upon our model. 

Detectable signals in the V spectrum (5 rms = 0.9 mJy) 
are restricted to the stronger of the two narrow peaks. 
This narrow peak is somewhat asymmetric and requires 
two Gaussians for a good fit. These two components 
(numbers 7 and 8) have field strengths of 10.9 ± 1.7 and 
17.9±0.9 mG, respectively. Three other Gaussians, num- 
bers 1, 6, and 11, have fields that are nearly 3er; however, 
visual examination of the V spectrum shows bumps and 
wiggles throughout at the « 1 mJy level, and we have 
no confidence in these purported fields. 

Since the significant splittings occur where the Stokes 
/ profile has the highest flux density, and since the 1667 
MHz transition dominates in all OHMs, we feel comfort- 
able assuming that the emission is from the 1667 MHz 
transition. Of course, because of the blending of the 1665 
and 1667 MHz lines, there is an unresolvable ambiguity 
that could affect the derived field strengths. 

The Stoke s I line profile has changed si nce the source's 
discovery by iDarling fc Giovanellil (|2001l ) . The flux den- 
sity of the broad component and the narrow component 
near 1372.3 MHz have remained the same. However, the 
narrow component near 1371.3 MHz, which used to be 
roughly 6 mJy (using our classical definition of Stokes 
/) weaker than their 32.5 mJy peak at 1372.3 MHz, has 
flared and is now the strongest component with a flux 
density of 44 mJy. This time-variable component is the 
same one that exhibits the largest splitting and there- 
fore probes the strongest field; in §[H]we compare this re- 
sult with newly-observed strong field detections in time- 
variable Galactic OH maser components. 

5.1.6. IRAS F12112+0305 

We were able to fit the Stokes / OHM emission quite 
nicely with 5 Gaussians as shown in Figure [6] The fit 
parameters are listed in Table [5] There is no detectable 
signal in Stokes V (S lms = 1.2 mJy). This is the first 
published spectrum of OHM emission in 12112. 

5.1.7. IRAS F14070+ 0525 

Table [7] lists the parameters for the 7 Gaussian com- 
ponents used to fit the Stokes / OHM emission in 14070. 
Since the 1665 and 1667 MHz lines are clearly blended 
in this source, we assume that all of the components rep- 
resent 1667 MHz emission. As seen in Figure [3 this 
decomposition provides a decent fit but there is no de- 
tectable signal in Stokes V {S nnB = 0.5 mJy). Gaussian 
number 1 shows a nearly 3<r detection of magnetic field; 
however, the associated feature in the Stokes V spectrum 
appears to be no more significant than the other features 
of mjy-strength intensity. We have no confidence in this 
near detection. 



5.1.8. IRAS F15327+2340 (Arp 220) 

The line profile of the Stokes I OHM emission in Arp 
220 is very complex. Our fit required 18 Gaussian com- 
ponents, as seen in Figure [51 to obtain reasonable resid- 
uals. An 18-component fit might seem overwhelming, 
but the components are easily obtained using our selec- 
tion guidelines from § 14.41 There are 9 distinct peaks 
(i.e., local maxima), represented in Figure[5]by Gaussian 
components 0, 4, 5, 7, 9, 13, 14, 16, and 17. There are 6 
narrow or fairly narrow bumps or shoulders that are not 
intense enough to produce local maxima, represented by 
Gaussian components 1, 2, 3, 11, 12, and 15. Compo- 
nent 10 was needed to represent the asymmetry in the 
brightest peak near 1638.15 MHz. Finally, components 
6 and 8 were needed to represent core-halo structure in 
the overall profile. This 18-component fit reproduces all 
of the visually obvious narrow, weak bumps, as well as 
the overall profile shape. However, the residuals exhibit 
a different signature in the line from that off the line, 
which means that our fit does not represent the I profile 
perfectly. We expended considerable effort making sure 
that each of the 18 Gaussians listed in Table[5]is actually 
needed for the fit by inspecting the residuals for different 
combinations of omitted Gaussians. 

Six narrow Gaussians exhibit visually obvious signa- 
tures in Stokes V (which has an rms noise of SVms = 1.16 
mJy) and provide good fits for Zeeman splitting. The ab- 
solute values of the derived field strengths range from 0.7 
to 4.7 mG, with 4 negative and 2 positive fields. Gaus- 
sians 2 and 3 have opposite field strengths. 

Gaussian numbers 9 and 11 are strong (a few hundred 
mJy), have comparable FWHMs of about 0.1 MHz, and 
are separated by about the FWHM. This makes them 
easily distinguishable. They have opposite field direc- 
tions as given by the least-squares fit, but the reversal 
in sign is also visually apparent. Zeeman splitting pro- 
duces a Stokes V pattern that looks like the frequency 
derivative of the line, with amplitude and sign scaled 
by Bt\. Thus, for a single Gaussian component, the V 
pattern looks like the letter "S" lying on its side, with 
an inevitable negative and positive part; the integral over 
Stokes V must be zero. However, the V pattern for these 
Gaussians in Figure [SJ doesn't look like this; instead it is 
positive on both sides of the line and negative in the mid- 
dle. The only way to obtain positive V on each side of a 
spectral bump i s for the field to have different signs on 
the two sides (cf . I Verschuurll 19691 Figure 3, the first radio 
detection of Zeeman splitting in Cas A). The integral of 
Stokes V must again be zero, with the central negative 
portion balanced by the two positive ones on the sides. 
The reversed field is not only a result of the fits, but is 
also visually apparent. 

We can compare our single-dish spectru m with the se- 
lected glob al VLBI spectr a pres ented by iRovilos et al.1 
( 2003) and lLonsdale et alj (| 19981) . There are a number 
of Gaussian components that appear to be directly asso- 
ciated with the resolved OHM spots: component 11 at 
5334 km s _1 originates in a southwestern spot tracing a 
positive field; component 6 at 5393 km s _1 originates in 
the southeast and traces a positive field; component 4 
at 5425 km s _1 originates in the center of the northeast 
ridge and traces a negative field; component at 5533 km 
s _1 originates in one of the southwestern spots tracing 
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a negative field. Three other features are more ambigu- 
ous: components 2 and 3 could be associated with either 
the northwestern or northeastern OHM ridges, while the 
brightest component, number 9 at 5351 km s _1 , appears 
to contain emission from both the northern ridges as well 
as the southwestern maser spots; these ambiguities pre- 
vent any possible field associations. The picture painted 
by the possible associations is for a field reversal from 
positive to negative from the southern to the northern 
features of the eastern OHM spots; there is no obvious 
reversal in the western region, but it is possible given the 
associations above. 

5.2. Linear Polarization 

For all observations, we used dual-polarized feeds with 
native linear polarization. This means that the observed 
Stokes U bs and V bs come from cross-correlation prod- 
ucts, which insulates them from system gain fluctuations. 
However, Stokes Qobs comes from the difference between 
the two native linear polarizations, so it is susceptible 
to time-variable, unpredictable gain fluctuations. This 
leads to coupling between Stokes / and Stokes Q b s ', in 
other words, a scaled replica of the I profile appears in 
the Qobs profile, with a random and unknown scaling 
factor, so Qobs is unreliable. 

Normally, when deriving linear polarization, one com- 
bines Stokes Qobs and U bs in the standard ways to obtain 
polarized intensity and position angle for the astronomi- 
cal source. However, since Q b s is unreliable for our mea- 
surements, we derived Stokes Q sr c and U src for the source 
from U bs alone by least-squares fitting its variation with 
parallactic angle. This is quite feasible at Arecibo be- 
cause all sources pass within 20° of the zenith, so track- 
ing for a reasonably long time provides a wide spread in 
parallactic angle. This makes the least-squares fit robust 
and provides good sensitivity and low systematics. For 
the source III Zw 35, which was plagued by serious inter- 
ference, we performed a minimum-absolute-residual-sum 
(MARS) fit. 

As with the Stokes I and V spectra, the least-squares 
derived Q src and U src spectra are displayed after sub- 
tracting both the off-source position and the continuum. 
We then use these baseline-subtracted Stokes spectra to 
derive spectra for polarized intensity and position angle. 
We do this because, even for the position-switched spec- 
tra, the continuum linear polarization is usually domi- 
nated by the diffuse Galactic synchrotron background. 
Although this prevents us from reliably deriving linear 
polarization for the (U)LIRG continuum radiation, the 
frequency-variable polarization is reliable. 

For two sources below, we least-squares fit for the Fara- 
day rotation measure RM. Performing this fit requires 
some care because the RM is derived from the position 
angle ip, which in turn is obtained by combining Q src 
and U src , which combine nonlinear ly through the arctan 
function [tp — 0.5 arctan(Q src /[/s rc )]. The channel-by- 
channel data are too noisy to produce a good-looking 
spectral plot of if), so on our plots we boxcar smooth by 
an appropriate number of points. One cannot linearly 
fit the unsmoothed values of ip to frequency because the 
arctan function produces nonlinear noise in ip. To avoid 
this problem, we performed a nonlinear fit to the un- 
smoothed arctan(Q src /[/ src ); this extra complication en- 
sures that the derived values and errors are unaffected by 



smoothing. 

We were unable to analyze the linear polarization for 
the two sources we observed using the GBT (IC 2545 and 
11506) because of inadequate parallactic angle coverage. 
We report the linear polarization for the Arecibo results 
here and discuss their interpretation in § 16.31 

We present three vertically stacked plots for each 
source below (Figures l9 Hl4|) . The top panel shows the 
position-differenced, baseline-subtracted Stokes / spec- 
trum over 12.5 MHz, therefore including both the 1665 
and 1667 MHz transitions for each source. The middle 
panel presents the linear polarization intensity and the 
bottom panel displays the derived position angle ip as a 
function of heliocentric frequency. 

5.2.1. IRAS F01417+1651 (III Zw 35) 

The linear polarization results for III Zw 35 are pre- 
sented in Figure [9l The top panels exhibit the position- 
differenced, baseline-subtracted Stokes / spectrum. The 
1665 MHz transition is clearly visible and the hypcrfinc 
line ratio is Rr = 6.0. The middle panels clearly display 
that the spectrum of linearly polarized intensity is ex- 
tremely spikey. Although higher S/N would help, the 
spikes might be real and possibly correspond to indi- 
vidual masers that are too weak to be seen clearly in 
the Stokes I spectrum. The polarized intensity shows a 
seemingly real peak centered near 1622.8 MHz, which is 
also the center of the Stokes I peak. The polarized in- 
tensity is about 5 mJy and the Stokes I peak is roughly 
500 mJy, so the fractional polarization w 1%. If the 
other spikes are real, then their fractional polarizations 
are much higher. 

The bottom panels display the position angle ip. Posi- 
tion angles exhibit less scatter than intensities, and the 
angle looks well-defined for the 1622.8 MHz peak. Also, 
it seems to show a gradual change across the line, which 
is about 1 MHz wide. The dashed line displays the result 
of a least-squares fit to the frequency dependence of the 
position angle, using only those points that are marked 
as diamonds: RM = -21900 ± 3700 rad m~ 2 . The ex- 
trapolated dashed line goes through the clusters of points 
associated with spikes centered near 1623.7 and 1624.0, 
and moreover, even the slope of the line matches the 
data for these spikes. The slope also seems to match the 
1624.4 MHz polarized-intensity spike, but the data are 
offset by about 60°. We speculate that: (1) these three 
polarized-intensity spikes come from individual OHMs 
that are too weak to see in the top panels of Figured (2) 
they all suffer the same Faraday rotation of ~ —21900 
rad m -2 as the central peak; and (3) the intrinsic po- 
sition angle for the 1624.4 MHz maser differs from the 
other two by about 60°. 

5.2.2. IRAS F10173+0829 

Figure \W\ displays the linear polarization results for 
10173. The polarized intensity shows a low-S/N spike 
that is centered on the Stokes / line: the polarization 
fraction is about 1% and the position angle about 60°. 
The spike is too narrow to fit for Faraday rotation. The 
hyperfme line ratio for 10173 is Rh = 10.7. 

5.2.3. IRAS F12032+ 1707 

The linear polarization results for 12032 are shown in 
Figure [Til The polarized intensity shows multiple spikes 
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that might be real. The most significant is centered at 
~ 1372 MHz, with a peak flux density of ~ 5 mJy, and 
has ip ~ 15°; near this frequency, I varies from ~ 10 
to ~ 40 mJy, so if this peak is real then the fractional 
polarization is huge, 50% to as 10% — unheard of for 
OH masers of any stripe. 

5.2.4. IRAS F12112+0305 

Figure rT2] shows the linear polarization results for 
12112. The lower frequencies are plagued by RFI, which 
remarkably disappears at the low-frequency boundary of 
the 1667 MHz line (centered at 1554.5 MHz). Accord- 
ing to the National Telecommunications and Information 
Administration Manual of Regulations and Procedures 
for Federal Radio Frequency Management, this RFI is 
likely attributable to space-to-Earth aeronautical mobile 
satellite communications operated by Inmarsat. There 
is no trace of any detectable linear polarization for this 
source. This is the first detection of the 1665 MHz tran- 
sition for 12112; the hyperfine line ratio is Rjj = 4.0. 

5.2.5. IRAS F 14-070+ 0525 

Figure [13] displays the linear polarization results for 
14070. The linear polarization intensity is approximately 
4 mJy across the entire 12.5 MHz bandwidth with an 
estimated position angle of —48°. 

5.2.6. IRAS F15327+2340 (Arp 220) 

The top panels of Figure [TJ] show the Stokes I profile 
for Arp 220 including both the 1665 and 1667 MHz tran- 
sitions. The hyperfine line ratio is Rh — 3.5. The middle 
panels show the linear polarization intensity, which has 
a well-defined peak centered at 1638 MHz and peaks at 
about 2 mJy. This is only 0.3% of the total intensity 
at this frequency. This is a very small fractional polar- 
ization but is very well-detected. 

The bottom panels show that the position angle of 
linear polarization is well-defined in two regions of low 
noise, one centered near 1638 MHz and the other near 
1636 MHz. The former region corresponds to the 1667 
MHz line and the latter is aligned with the 1665 MHz 
transition. This 1665 MHz line is unconvincingly visible 
in the polarized intensity spectrum, but the low noise in 
its position angle spectrum is unmistakable. 

We fit the frequency variation of tp to obtain the Fara- 
day rotation measure RM using those points marked 
as diamonds in the bottom right panel of Figure [Til 
For the 1638 MHz component alone, we obtain RM = 
5230 ± 7930 rad m~ 2 . For the combination of the 1636 
and 1638 MHz components, we obtain RM = 1250±1040 
rad m~ 2 . These errors are considerable and make the for- 
mal result only marginally significant. The dashed line 
in the bottom-right panel displays the result of the fit 
for both components together; visual inspection shows 
that not only is it an acceptable fit for both components 
together, but it is also acceptable for the 1638 MHz com- 
ponent alone. It is not unreasonable to conclude that the 
OHM radiation from both OH lines suffers a common 
Faraday rotation of RM w 1250 rad m~ 2 ; this is ~ 20 
times smaller than the value derived for III Zw 35. 

6. DISCUSSION 
6.1. OH Maser Zeeman Pairs in the Milky Way 



While our results are the very first in situ Zeeman 
detections in external galaxies, OH masers in the MW 
have been used as Zeeman magnetometers for well over 
a decade. In contrast to OHMs, Galactic OH maser emis- 
sion lines are so narrow (~ 0.5 km s" 1 ) that fields of « 1 
mG are sufficient to completely split the left and right 
circular a components into pairs. Mor e than 100 of these 
Zee man pairs have been com piled by iFish et all (|2003l ) 
and lReid &; Silversteinl (|1990f ) with a distribution whose 
mean is consistent with mG and whose standard de- 
viation is 3.31±0.09 mG. Typical densities in OH maser 
regions are n ~ 10 6 — 10 7 cm -3 ; for a field strength of 
~ 10 (iG in gas at ~ 1 — 100 cm" 3 , the fields probed by 
Galactic OH masers are consist ent with the enhancement 
of \B\ oc n 1 / 2 (|Fish et al.ll2003[ ). The linear polarization 
of the a components is often measured in addition to the 
7r component, but the it components, which are in the- 
ory 100% linearly polarized, are rarely measured to be 
purely so. 

Unlike in the OH masers in our Galaxy, the flux density 
of the 1667 MHz transition in all OHMs is larger than 
that of the 1665 MHz transitio n and , until now, no polar- 
ization has been detected (|Ldl2005l ). There is no defini- 
tive explanation for the dominance of the 1667 MHz tran- 
sition, but recent work suggests that this probably arises 
because the extragalactic lines are wider than the Galac- 
tic maser lines (P. Goldreich, privat e communication; 
M. E litzur, private communication; lLockett fc E litzur 

Hool . 

Our detections yield a median line-of-sight magnetic 
field strength of ~ 3 mG in OHMs in (U)LIRGs, which is 
comparable to the field strengths measured in OH masing 
regions in the MW. This strongly suggests that the local 
process of massive star formation occurs under similar 
conditions in (U)LIRGs, galaxies with vastly different 
large-scale environments than our own. 

The magnetic field strengths we find in the OHMs 
in (U)LIRGs (~ 3 mG) are comparable to the volume- 
averaged fields of > 1 mG inferred from synchrotron ob- 
servations. These results imply that mG magnetic fields 
likely pervade most phases of the interstellar medium in 
(U)LIRGs. It is unclear, however, how to physically re- 
late the two different magnetic field strengths in more 
detail given the possibility that each may probe rather 
different phases of the ISM. Some models of OHMs invoke 
radiative pumping in mole cular clouds with gas densities 
- 10 3 - 5 - 10 4 cm- 3 (e.g.. iRandell et al.lll995l ). This is 
similar to the mean gas density in the central ~ 100 pc 
in (U)LIRGs, in which case our observations likely probe 
the mean ISM magnetic field (whether the synchroton ra- 
diation also arises from gas at this density is unclear; up- 
coming GLAST ob servations of neutral p ion decay may 
help assess this; see lThompson et aLl 2007). It is also pos- 
sible, however, that the OHMs arise in somewhat de nser 
gas (n ~ 10 6 — 10 7 cm" 3 ; e.g.. lLonsdale et aT, ll998h . as 
appears to be true in the MW fe.g.. lFish et al. [20031 ). In 
this case, the magnetic field probed by OHMs is likely 
stronger than that in the bulk of the ISM. If we as- 
sume the B oc n 1 ^ 2 scaling often a ssumed in the MW 
(|Mouschoviadll976t [Fish et al.ll2003f) . the field strengths 
in the masing regions in (U)LIRGs are probably within 
a factor of ~ 3 of the mean ISM field (rather than a 
factor of several 100 in the MW), given the large mean 
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gas densities in (U)LIRGs. This is still reasonably con- 
sistent with the mean field strength of > 1 mG inferred 
from synchrotron observations. Without a better under- 
standing of the physical conditions in the masing regions, 
however, it is difficult to provide a more quantitative 
connection between our inferred field strengths and ei- 
ther the mean ISM field or the magnetic field probed by 
synchrotron emission. Ultimately doing so is important 
because it will allow stringent constraints to be placed 
on the dynamical importance of magnetic fields across a 
w ide range o f physi cal conditions in (U)LIRGs. 

iFish et ail ((2003), with their comprehensive survey 
of Galactic OH masers and the accompanying statisti- 
cal discussion, strongly support several previous sugges- 
tions that the field direction in OH masers usually mir- 
rors that of the large-scale field in the vicinity of the 
m asers. MERLIN ob s ervat ions of OH masers in Cep A 
by iBartkiewicz et~a l. ( 2005) also present Zeeman detec- 
tions corroborating the field's alignment with the ambi- 
ent ISM field direction. Thus, measuring the direction 
of the field in an OH maser reveals the field direction 
not only in the maser, but also outside and in the vicin- 
ity of the OH maser. For the MW, this aids us to infer 
the large-scale magnetic field morphology. To directly 
compare the star formation processes in the MW and 
(U)LIRGs, it will be necessary to increase the sample 
of magnetic field strengths in (U)LIRGs and to directly 
map the Zeeman splitting of individual OHM spots us- 
ing VLBI in order to probe whether reversals occur at 
smaller angular scales. 

6.2. Strong Fields and Time Variability 

ISlvsh fe Migenei ([2006] ) and lFish fe Reidl (f2007h both 
observed fields of 40 mG using Zeeman observations of 
OH maser spots in W75N; these are the highest field 
strengths measured in Galactic OH masers and are an 
order of magnitude larger than the typical OH maser 
field. These OH maser spots also happen to have been 
flaring based on multi-epoch VLBA observations; per- 
haps time variability in OH masers is correlated with 
strong magnetic fields. Interestingly, our strongest de- 
tection, Bn ~ 18 mG in the gigamaser 12032, occurs in 
an OHM component that has increased in flux density by 
a factor of two since its pr evious published observation 
([Darling fe Giovanellill200l . 

These results strongly support the development of an 
observational program to monitor both the time vari- 
ability of the Stokes / flux density and magnetic field 
strength in OHMs as well as the necessity of observing 
the circular polarization of time- variable Galactic masers 
in hopes of detecting strong magnetic fields. 

6.3. Linear Polarization and Faraday Rotation 

Our measured rotation measures of RM ~ 21900 rad 
m -2 for III Zw 35 and RM ~ 1250 rad m~ 2 for Arp 
220 are large by most standards, but are not unreason- 
able for (U)LIRGs. As mentioned in § [TJ the magnetic 
field strength throughout the ULIRG ISM should be > 1 
mG from synchrotron observations. Electron densities 
are estimated to be ~ 1 — 10 cm~ 3 in the hot ion- 
ized plasma, both from observations of X-ray emission 
(e.g. jGrimes et al"1l2005f ) and from th eoretical models of 
supernova-driven galactic winds fe.g.. lChevalier fe Cleggj 



1985). Over a path length ~ 100 pc in the central por- 
tions of ULIRGs, n e ~ 1 cm -3 and B ~ 1 mG imply 
(n e B\\L) - 0.1 G cm- 3 pc, or RM ~ 80000 radians 
m~ 2 . This is a factor of 4-60 larger than our measured 
values. 

It is reasonable for this simple estimate to overestimate 
the measured RM. This is because the RM depends 
only on the line-of-sight field component. The proba- 
bility density function for the line-of-sight component of 
a randomly oriented magnetic field is flat between zero 
and the perfectly-oriented case; thus, for a set of sources 
with randomly-oriented fields, the observed line-of-sight 
field component is reduced by a factor of two, and 1/4 
of the sources have the observed component less than 
1/4 the perfectly-aligned value. In addition, and prob- 
ably more importantly, the observed Faraday rotation 
responds only to the systematic line-of-sight field com- 
ponent, while the synchrotron radiation and the total 
magnetic energy depend on the total field, systematic 
plus random. Our estimate of RM ~ 80000 is for the 
total field, not the systematic field, because the latter is 
much harder to predict. 

The measured RM might also be reduced by finite 
source-size effects a nd / or prop agation through an inho- 
mogenous medium (|Burn|[l966T ). First, suppose that the 
magnetic field is everywhere uniform but that the Fara- 
day rotation is produced in the same region where the 
maser radiation is produced, and that this region is ex- 
tended along the line of sight. In this case, different 
line-of-sight depths of the maser are rotated by different 
amounts. This washes out the linear polarization and 
can reduce the apparent Faraday rotation. In the other 
extreme, think of the field as primarily random except for 
a small uniform component. Maser radiation observed at 
a given frequency might come from more than one maser 
located at different positions on the sky or at different 
distances into the source. In the former case, the RM 
might change with position on the sky; in the latter, it 
might change along the line of sight. In cither case, its 
average value can be small. In addition, for an individ- 
ual maser the field might fluctuate along the line of sight, 
reducing the total RM. 

The interpretation of the linear polarization and RM 
data is thus currently difficult and non-unique. Obser- 
vations of more systems would be helpful and may ulti- 
mately provide unique constraints on the thermal elec- 
tron density and/or magnetic field structure (e.g., rever- 
sals) in the nuclei of (U)LIRGs. 
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TABLE 1 

IRAS F01417+1651 (III Zw 35) Gaussian Fit Parameters 
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Gaussian 


(mJy) 


(MHz) 


(MHz) 


(km s _1 ) 


(mG) 


(1) 
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(3) 


(4) 


(5) 


(6) 





190.71 ±6.77 


1622.3743 ±0.0049 


0.2041 ±0.0081 


8312.6 


2.94 ±0.18 


1 


255.68 ±22.16 
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8253.0 


-0.47 ±0.18 
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99.26 ±6.28 
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1.73 ±0.78 


3 


176.13 ±34.75 


1622. 7604 ±0.0020 


0.0905 ± 0.0080 


8239.3 


-2.73 ±0.13 


4 


1.17 ±59.50 


1622.8864 ±0.0142 


0.1544 ± 0.0239 
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-3.59 ±0.26 



Heliocentric Velocity (km s ') 
8400 8300 8200 8100 




1621.5 1622.0 1622.5 1623.0 1623.5 1624.0 

123 4 



10 



I 

> 



-10 





i | i i i i | i i i i | i i i i 






y - 









1621.5 



1622.0 1622.5 1623.0 1623.5 

Heliocentric Frequency (MHz) 



1624.0 



Heliocentric Velocity (km s ') 
8400 8300 8200 8100 



500 F" 
400 r 
g 300 - 

s 

m 200 ' 7 
100 - 




1621.5 



1622.0 



1622.5 



1623.0 



1623.5 




1624.0 



1622.0 



1622.5 



1623.0 



1623.5 




1624.0 



1621.5 



1622.0 1622.5 1623.0 1623.5 

Heliocentric Frequency (km s ') 



1624.0 



Fig. 1. — Total intensity and circular polarization results for IRAS F01417+1651 (III Zw 35). (Top left) The Stokes / spectrum (solid 
line; twice the conventionally defined flux density) is plotted as a function of heliocentric frequency (bottom axis) and optical heliocentric 
velocity (top axis). The profile of each Gaussian component is plotted as a dashed line with the corresponding component number shown 
below the frequency axis at the corresponding central frequency. Residuals from the composite Gaussian fit are plotted through the center 
of the panel (solid, light-weight line) and are expanded by a factor of 4. The scale bars near the right edge of the plot correspond to the 
labelled flux density range. (Bottom left) The Stokes V spectrum (solid line) and its fit (dashed line). (Top right) Composite Gaussian fit 
to Stokes /. (Middle right) Measured Stokes V (solid line); hypothetical Stokes V (dashed line) produced by a uniform By = 1 mG using 
the derivative of the composite / profile above. (Bottom right) Derived By (crosses and solid line) and uncertainty (error bars) from the 
B(v) fit. All spectra and residuals are smoothed with a boxcar of 7 channels. 
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TABLE 2 

IRAS F10038-3338 (IC 2545) Gaussian Fit Parameters 



S v Av v & B\\ 

Gaussian (mJy) (MHz) (MHz) (km s" 1 ) (mG) 

(1) (2) (3) (4) (5) (6) 



47.06 ±3.10 1612.3870 ± 0.0044 0.1505 ± 0.0125 10221.0 1.58 ±0.75 

1 179.73 ± 5.15 1612.8661 ± 0.0032 0.1747 ± 0.0077 10128.9 - 1.76 ±0.26 

2 83.85± 5.50 1613.0090± 0.0119 0.9548 ± 0.0405 10101.4 -11.28 ± 1.16 

3 227.46 ± 36.43 1613.0362 ± 0.0052 0.1222 ± 0.0095 10096.2 - 0.07 ±0.18 

4 2.85 ± 14.04 1613. 1570 ± 0.0089 0.1703 ± 0.0157 10073.0 - 0.13 ±0.16 

5 97.85 ±4.75 1613.3556 ± 0.0035 0.1136 ± 0.0075 10034.9 1.67 ±0.33 
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Fig. 2. — Total intensity and circular polarization results for IRAS F10038— 3338 (IC 2545). See caption for Figure Q] (Top) Residuals 
are expanded by a factor of 4. 
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TABLE 3 

IRAS F101 73+0829 Gaussian Fit Parameters 



S v Av v Q B\\ 

Gaussian (mJy) (MHz) (MHz) (km s" 1 ) (mG) 

(1) (2) (3) (4) (5) (6) 



30.62 ± 6.34 1589.2420 ± 0.0037 0.0518 ± 0.0113 14735.9 2.34 ±2.50 

1 16.98 ± 5.95 1589.2862 ±0.0041 0.0244 ± 0.0106 14727.1 4.19 ±3.42 

2 161.91 ± 7.80 1589.3190 ± 0.0024 0.1786 ± 0.0050 14720.7 0.25 ±0.89 

3 15.54 ± 7.63 1589.3258 ± 0.0073 0.0463 ± 0.0240 14719.3 -2.95 ±5.15 

4 0.42 ± 1.58 1589.5315 ±0.0194 0.5661 ± 0.0257 14678.6 0.93 ±5.54 

5 3.50 ± 2.82 1589.6383 ±0.0154 0.0411 ± 0.0412 14657.5 -2.90 ±17.78 

6 16.62 ± 2.46 1589.8772 ± 0.0072 0.1670 ± 0.0279 14610.2 0.80 ±7.57 
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Fig. 3. — Total intensity and circular polarization results for IRAS F10173+0829. See caption for Figure [T] (Top left) Residuals are 
expanded by a factor of 5. (Middle right) Dashed line shows the hypothetical Stokes V produced by a uniform By = 2 mG using the 
derivative of the composite / profile above. 
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TABLE 4 

IRAS Fl 1506-3851 Gaussian Fit Parameters 



S v Av vq £?|| 

Gaussian (mJy) (MHz) (MHz) (km s" 1 ) (mG) 

(1) (2) (3) (4) (5) (6) 



64.73 ± 20.75 1650.0094 ± 0.0428 0.1910 ± 0.0456 3152.3 1.21 ±0.27 

1 102.55 ± 26.92 1650.1891 ± 0.0240 0.1810 ± 0.1018 3119.3 0.45 ± 0.20 

2 86.37 ± 83.80 1650.3232 ± 0.0231 0.1358 ± 0.0467 3094.7 0.36 ±0.17 

3 116.30 ± 24.97 1650.4685 ± 0.0257 0.2693 ± 0.0976 3068.0 0.73 ±0.20 

4 0.31 ± 7.25 1650.7843 ±0.1171 0.3616 ± 0.2594 3010.1 0.68 ±0.75 

5 23.79 ± 7.46 1651.0731 ± 0.0102 0.1260 ± 0.0368 2957.1 1.03 ±0.50 
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Fig. 4. — Total intensity and circular polarization results for IRAS F11506— 3851. See caption for Figure[T] (Top) Residuals are expanded 
by a factor of 2. 
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TABLE 5 

IRAS F12032+1707 Gaussian Fit Parameters 
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Fig. 5. — Total intensity and circular polarization results for IRAS F12032+1707. See caption for Figure[T] (Top) Residuals are expanded 
by a factor of 2. 
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TABLE 6 

IRAS F121 12+0305 Gaussian Fit Parameters 



S v Av v & B\\ 

Gaussian (mJy) (MHz) (MHz) (km s" 1 ) (mG) 

(1) (2) (3) (4) (5) (6) 



16.19 ± 0.75 1554.1822 ±0.0013 0.1467 ± 0.0054 21831.1 -0.20 ±1.99 

1 58.16 ±0.51 1554.2767 ±0.0048 0.5188 ± 0.0098 21811.6 0.25 ±1.17 

2 48.63 ± 0.97 1554.5634 ± 0.0007 0.1893 ± 0.0028 21752.3 1.02 ±0.77 

3 13.83 ± 0.59 1554.5641 ±0.0006 0.0309 ± 0.0017 21752.1 0.27 ±1.03 

4 0.32± 0.34 1554.8544 ±0.0037 0.4359 ± 0.0063 21692.1 3.09 ±1.92 
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Fig. 6. — Total intensity and circular polarization results for IRAS F12112+0305. See caption for Figure[T] (Top) Residuals are expanded 
by a factor of 4. 
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TABLE 7 

IRAS F14070+0525 Gaussian Fit Parameters 



S v Av vq £?|| 

Gaussian (mJy) (MHz) (MHz) (km s" 1 ) (mG) 

(1) (2) (3) (4) (5) (6) 



3.50 ±0.22 1315.2359 ± 0.0105 0.4088 ± 0.0317 80262.3 -10.42 ±10.49 

1 8.95 ±0.19 1316.0759 ± 0.0715 2.1625 ± 0.0861 80019.7 26.65 ±9.57 

2 5.70 ±0.30 1316.2886 ±0.0047 0.2040 ± 0.0134 79958.4 1.76 ±4.65 

3 11.26 ± 0.68 1316.9061 ±0.0100 0.6771 ± 0.0302 79780.3 - 3.31 ±4.67 

4 0.08 ±0.50 1317.5519 ±0.0118 0.4878 ± 0.0271 79594.2 - 3.16 ±5.59 

5 2.78 ±0.26 1318.8802 ± 0.0106 0.2664 ± 0.0307 79212.1 12.72 ±10.82 

6 10.76 ±0.16 1319.1605 ±0.0113 1.2737 ± 0.0199 79131.6 -12.92 ±6.16 
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Fig. 7. — Total intensity and circular polarization results for IRAS F14070+0525. See caption for Figure[T] (Top) Residuals are expanded 
by a factor of 2. 



Extragalactic Zeeman Detections in OH Megamasers 



TABLE 8 

IRAS F15327+2340 (Arp 220) Gaussian Fit Parameters 
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Fig. 8. — Total intensity and circular polarization results for IRAS F15327+2340 (Arp 220). See caption for Figure[T] (Top) Residuals 
are expanded by a factor of 16. 
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Fig. 9. — Linear polarization results for IRAS F01417+1651 (III Zw 35). (Left) Top panel shows Stokes /, middle panel shows the linearly 
polarized intensity, and bottom panel shows the position angle over the entire 12.5 MHz bandwidth. (Right) Same as left panels, but with 
the frequency range narrowed to 6 MHz; the bottom right plot also shows the fitted Faraday rotation as a dashed line whose slope was 
determined by fitting to the points marked as diamonds. All spectra are plotted as a function of heliocentric frequency (bottom axis). The 
top panels show the optical heliocentric velocity (top axis). All spectra are smoothed by a boxcar of 23 channels. 
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Fig. 10. — Linear polarization results for IRAS F10173+0829. See caption for Figure [9] All spectra are smoothed by a boxcar of 17 
channels. 
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Fig. 11. — Linear polarization results for IRAS F12032+1707. See caption for Figure [9] All spectra are smoothed by a boxcar of 31 
channels. 
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Fig. 12. — Linear polarization results for IRAS F12112+0305. See caption for Figure [9] All spectra are smoothed by a boxcar of 11 
channels. 
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Fig. 13. — Linear polarization results for IRAS F14070+0525. See caption for Figure [9] All spectra are smoothed by a boxcar of 11 
channels. 
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Fig. 14. — Linear polarization results for IRAS F15327+2340 (Arp 220). See caption for Figure[9] (Right) The frequency range has been 
narrowed to 4 MHz. All spectra are smoothed by a boxcar of 9 channels. 



